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The Pais-Uhlenbeck(PU) oscillator is the simplest model with higher time derivatives. Its prop¬ 
erties were studied for a long time. In this paper, we extend the 4th order free PU oscillator to a 
more non-trivial case, dubbed the 4th order time dependent PU oscillator, which has time depen¬ 
dent frequencies. We show that this model cannot be decomposed into two harmonic oscillators in 
contrast to the original PU oscillator. An interaction is added by the coordinate transformation of 
Smilga. 


I. INTRODUCTION 

It is worth investigating higher derivative the¬ 
ories. A higher derivative theory usually vi¬ 
olates unitarity, and the Pais-Uhlenbeck(PU) 
oscillator [Tj is the easiest model for studies 
what the higher derivatives cause. In higher 
derivative theories, their Hamiltonian needed 
for quantization and examining stability can 
be obtained by Ostrogradski’s method[2, [lOl| . 
The Hamiltonian of the PU oscillator looks lit¬ 
tle complicated, but it can be rewritten to the 
very simple form of several independent har¬ 
monic oscillators. For example, as regards the 
4th order PU oscillator, whose equation of mo- 
tion(EOM) has up to the 4th order time deriva¬ 
tive, its Hamiltonian iJpu can be separated into 
two independent harmonic oscillator Hamilto¬ 
nians, i?pu = Hi — H 2 , by Smilga’s canoni¬ 
cal transformation^. In the last form, quanti¬ 
zation is quite simple, just by quantizing each 
of harmonic oscillators. It is obvious that the 
negative sign before H 2 means that the en¬ 
ergy is not bounded from below, so that this 
model may be unstable (once interactions are 
added) d, The negative energy generically 
comes out in higher derivative theories Q, how¬ 
ever, for the PU oscillator model no instabil¬ 
ity occurs because it is a free modelfi, 0, 0-111 ■ 
Given some interactions, maintaining separabil¬ 


ity is known to play the key role for the stability 
of the interacting PU model. It is significant to 
consider other interacting higher derivative PU 
models, and examine whether the separability is 
conserved or not. In this paper, we extend the 
4th order PU oscillator to the time dependent 
4th order PU(tdPU) oscillator, whose frequen¬ 
cies depend on time, to see whether it has the 
separability or not. Our model may be consid¬ 
ered as a PU oscillator under the influence of 
some external forces, so that the energy is not 
conserved. In this situation, when one naively 
uses Smilga’s transformation, a problem arises. 
The tdPU Hamiltonian obtained by Ostrograd¬ 
ski’s method and the Hamiltonian transformed 
by Smilga are in one to one correspondence with 
each other, however, their equations of motion 
are not. This means that Smilga’s transforma¬ 
tion is not canonical in our model, so that it is 
necessary to add some correction terms to the 
naively transformed Hamiltonian. Here we find 
those correction terms from the two points of 
view: (i) a comparison of the differential equa¬ 
tions, and (ii) the generating function method. 
We show that the separability of our model is 
deformed because the corrections become the 
interaction terms of two time dependent har¬ 
monic oscillators. 
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II. PREVIOUS RESULTS 

A. 4th order PU oscillator 

In EH , the 4 th order PU oscillator is defined 
as 

Lpu(a;, i, x) = (1) 

where Oi := ujiOJ2,^2 '■= ■ Let us 

introduce the PU Hamiltonian by applying Os- 
trogradski’s method^ to (P). The PU Hamil¬ 
tonian iLpu is 

Hpu{qi,q2,Pi,P2) = Piqi +P 292 - Lpu 

1 2 

= Pi 92 + 2^2 

+ ^^ 292 ^ - ^niqj^,{2) 
where qi,Pi (i = 1,2) are canonical coordinates, 

'qi = X, 

dLpu d dLpu 

dx dt die 

< = -fl2X — X , ( 3 ) 

(?2 = X, 

dLpTj 


to 

Lfpu(Qi,Q2,Pi,P2) =iLi(Qi,Pi) 

-iL2(Q2,P2), (5) 

where H^{Qi,Pi) := -I- 

Thus we see that the PU oscillator can be 
separated into two harmonic oscillators. Smilga 
used this property to quantize the PU oscillator 
and and investigate its quantum stability. 

III. OUR RESULTS 

A. Time dependent PU oscillator 

Our Lagrangian of the time dependent 
PU(tdPU) oscillator is 

LtdPu(x,x,x) = ix^ - U2(t)x^ 

-bUi(t)x^ (6) 

where Hi(t) := a;i^(t)w2^(f), ^2(0 := + 

uj2^{t). This Lagrangian can be obtained by re¬ 
placing uji —>• uji{t){i = 1 , 2 ) in (HJ. Its EOM 
is 

x'-f f22(t)x-I-H 2 (t)x-b Hi(t)x = 0. (7) 

The Hamiltonian and its EOM are 


B. Decomposition of the PU Hamiltonian 
by Smilga’s transformation 

In 0 , Smilga gave the transformation that 
separates Hp\j as follows: 


9 1 — —(wi (52 — Pi), 

Wi 

Pi = wi7(wiP 2 - a;2^Qi), 

92 = 7 (wiQi - P2), 

P2 = 7(wiPl - W2^Q2)- 


( 4 ) 


L 7 tdPu( 9 i, 92 ,Pi,P 2 ) = Pi 9 i +P292 — LtdPU 


= Pi 92 + 2^2 


+ 2 ^ 2 (i) 92 ^ 



= A 



( 8 ) 

( 9 ) 


where 7 := tj— UJ2 ■ Iii this paper, we 
call it Smilga transformation. 

Here we assume uji > UJ2 for simplicity, and 
this transformation is a canonical transforma¬ 
tion. According to (| 3 ]), Hpu can be rewritten 


where 

/ 0 1 0 0\ 

A- ^ 0 01 

fliit) 0 0 0 

V 0 - 02 (t) - 10 / 


( 10 ) 
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Here Pi,qi {i = 1,2) are defined by ([ 3 ]) with time 
dependent frequencies, 0 and 0 are equiva¬ 
lent in our case. 

By using Q with time dependent frequencies, 
([8]) can be written as 

^^tdPu(Qi, Q2, -Pi, ^2) 

= (Ptd)i(Qi,Pi) - {HMQ2.P2), ( 11 ) 

where (Ptd)i(Qi,Pi) = ^P*^ + ■ Its 

EOM is 


where 


^ ■" 0 0 0’ 

V 0 0 0/ 

It is straightforward to check that ([8]) is equal 
to (|TT 1 ) under the Smilga transformation. How¬ 
ever, the time dependent Smilga transforma¬ 
tion is not a canonical transformation. In other 
words, it doesn’t reproduce the EOM © or, 
equivalently, does not follow from (1121) . In 
the next two Sections, we show that some cor¬ 
rection terms are needed to correct the situa¬ 
tion, and they become interaction terms of the 
harmonic oscillators. 



( 12 ) 


1 0 
0 -1 




B. Comparison of EOMs 

Here we find the correction term by a compar¬ 
ison of the EOMs ([ 9 l) and (fT^ . Let us introduce 
some notation first. 

X := {qi,q 2 ,Pi,P 2 f, ( 14 ) 

r:=(Qi,Q2,Pi,P2)^, ( 15 ) 

Here ^ denotes transpose of a matrix. Y, X and 
A, B are related by 


where M is the (inverse) coefficient matrix of 
the Smilga transformation (| 4 ]). With these no¬ 
tation, the EOM m and m can be written 
as 

tA = 

^Y = BY + a. (19) 

Here a represents a correction to (fT^ . Differ¬ 
entiating the first relation (ITSl) with respect to 
t and substituting dUD to the left hand side, we 
hnd 

Ay = M—X + MX = MAX + MX 
at at 

= BY + a = BMX + a, 
a = MX = MM-^Y. 


Therefore, a ^ 0 , unless the model is time in¬ 
dependent. More explicitly, a is 


a = 'y^it) 


/Cl 0 0 El 

0 02 Pi 0 

0 P2 C3 0 
\ P2 0 0 C4 




where 
Cl := 
C2 := 
C 3 := 
C 4 := 


7 (i) , wi(t) 


73 (t) 
ijt) 
73 (t) 
i(t) 
73(f) 
i(.t) 

73(f) 


112(f), 


Wi(f) 

-|- 2uji (f)wi (f) 
Ci(f) 

wi(f) 72 (f) 

- 2w2(f)w2(f) 


— 2a;2(f)w2(f), 


El — 2 uji{t), P2 :— 2a;i(f)w2(f)w2(f), 
7(f) := l/ijuji^{t)-u}2^{t). 


Y = MX, 
MA = BM, 


(16) 

(17) 
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C. Generating function method 

In this part, we show that the correction 
term can be also obtained by using a canonical 
Hamiltonian The latter can be found by 

the generating function method, via solving the 
equations 

_ dW{q,Q,t) _ dW{q,Q,t) 

Pi r\ 5 r\ 5 

UQi 

TJ I ^ TJ! _ 9 W{q,Q,t) 

HtdV\j{(l^p) — ^tdPUvQ’ 7 

where q := {qi,q2),P ■= iPi,P2),Q ■= 

{Qi,Q 2),P ■= {Pi,P2) and W{q,Q,t) is the 
generating function of a canonical transforma¬ 
tion. Then we find 

W{q,Q,t) = -uj^{t)qiq2 + 

lit) 

H— 7 -rQ 292 ~ ‘^iit)QiQ2 + git)- 
lit) 

Here g{t) is an arbitrary time dependent func¬ 
tion. Using the above we find that the 
^tdPu(Q. P) is given by 

, dW 

PtdFviQ^P) — HtdPviQj P) + 

= HtAMQ,P) + Hinu (20) 

where 

Hint ■■= - 27 ^(t)a;i(<)w 2 (<)w 2 (t)QiQ 2 
-2-f\t)uJiit)PiP2 

+ ^ ^ i‘^lit)^2it) + Wi(t)w 2 (t))Pl( 5 l 

Wl(t) 

-b7^ (f ) (wi (t) Wi (t) -b W2 (<) W2 (t) ) P 2 Q 2 
+git) (21) 


where 

Fl = {uJl{t)uJl{t) + UJ2it)uJlit)), 

F 2 = UJi{t)uJi{t) + UJ2it)uJ2it), 

F3 = {uJl(t)L02it) +UJ2it)uJlit)), 

F4 = -(wi(t)wi(t) -b W2(t)w2(t))- 

It is easily confirmed that Ci = Fi, {i = 
1 , 2 , 3 , 4 ). This means that the new terms in 
(l 2 nil reproduce a. Hence, the interacting Hamil¬ 
tonian (uni) is canonical and the correction a of 
the EOM (IT^ is equivalent to the interaction 
terms in the canonical Hamiltonian 

IV. CONCLUSION 

In this paper, we extended the free 4 th order 
PU oscillator to the non-trivial case, the 4 th or¬ 
der tdPU oscillator, whose frequencies depend 
on time. We found that our model cannot be 
written down in the form of two harmonic os¬ 
cillators. This is because the Smilga transfor¬ 
mation is not canonical in our extended model, 
so that it was necessary to add some correc¬ 
tion terms to make it canonical. We obtained 
those corrections from the two points of view: 
(i) a comparison of the differential equations, 
and (ii) the generating function method. We 
showed that the correction terms in the tdPU 
EOM can be written as interaction terms in the 
canonical tdPU Hamiltonian. 
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Its EOM is 


Ar^BY + j^(t) 


( Fl 0 
0 F2 
0 E2 
\E2 0 


0 El 

El 0 
F3 0 
0 F4 
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